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Introduction 


This  is  a  discussion  of  methods  for  describing,  mathematically, 
flows  between  compartments  in  a  multi -compartment  system.  We  will  give 
the  conventional  theory,  based  on  the  solution  of  a  system  of  linear 
differential  equations;  we  will  also  give  a  theory  based  on  probability, 
viewing  the  system  as  a  collection  of  "states"  with  a  particle  moving 
from  state  to  state  with  certain  probabilities,  remaining  in  each  state 
a  random  time  with  an  exponential  distribution.  finally,  we  will  take 
still  another  approach,  again  based  on  probability  theory,  in  which  we 
consider  the  sojourn  time  of  a  particle,  that  is,  the  time  it  spends 
after  leaving  a  given  compartment  before  returning  to  that  compartment. 
This  last  approach  seems  to  be  potentially  the  most  useful,  as  it  fits 
well  the  practical  situation  in  which  an  experimenter  is  able  to  observe 
the  disappearance  of  tracer  amounts  of  a  substance  from  only  one  of 
several  biological  "compartments" .  He  is  then  interested  in  finding 
models  consistent  with  his  observed  disappearance  curve.  In  our  parlance, 
we  will  look  for  models  with  a  given  sojourn  time.  The  procedures  for 
finding  a  class  of  models  with  a  given  sojourn  time  from,  say,  the  first 
compartment,  will  be  easier  than  those  based  on  solving  the  system  of 
linear  differential  equations,  in  addition,  the  methods  will  be  more 
general,  and  apply  to  non-linear  situations  as  well. 

2.  Multi -compartment  systems  viewed  as  a  set  of  simultaneous  linear 
differential  equations 


It  is  customary  to  express  the  amount  of  tracer  in  each  of  the 
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compartments  of  a  multi -compartment  system  as  the  solution  to  a  system 
of  linear  differential  equations.  We  give  a  brief  outline  of  a  method 
for  solving  such  systems,  using  matrix  theory.  Define  the  vector 
t,  =  (x^t)  ,x2(t)  ....  ,xn(t))  .  Wo  want  to  solve 

f,  =  a,  s(o)  =  *0, 

where  A  is  an  n  *  n  matrix  of  constants,  and  is  diagonable;  the  dot 

means  differentiation  with  respect  to  t.  We  express  A  in  terms  of 

its  principal  idempotents,  A  =  r.  E.  +•••+  r  E  ,  with  r r  the 

11  m  m  (A  -  rh)  ™ 

distinct  characteristic  roots  of  A  and  E,  =  TT  - ^ — •  Then 


rnt  r„t  r  t 

t.'A  i  .  c  r  .  .  m  ... 

e  =  e  E.  +  e  E„  +• • •+  e  E 
12  m 


l  r.  -  r . 

jfi  i  J 


and  the  solution  to  (1)  is 


r1 1  r  t  r  t 

<  "  'o<e  E1  +  «  E?  +-+  '  "  £„)• 


As  examples,  consider  m  =  2  and  3*  When  A  has  two  distinct  roots, 
r^,Tp,  we  have 

A  -  r„I  r  t  A  -  r  I  r  t 

C  =  Sn[ - —  e  1  +  - i-  e  2  ]. 

0  rl  "  1 2  r2  -  rl 

If  A  has  distinct  roots  r^,r^,r^,  then  the  solution  to  (l)  is 

(A  -  r^l)  (A  -  r^l)  r-^t  (A  -  r^l)  (A  -  r^l)  r,,t 
^  ~  ^°^rl  "  rp)  (ri  "  r3^  '  +  (r2  -  ^(r.,  -  r^)  e 

(A  -  r1l)(A  -  r2l)  r^t 


r3  "  riMr3  -  r2; 


+ 


4 


or 


r  /2  4\  “6t  i  /I  4\  -2t 
t  =  (j.  -  |)e  +  (j,  |)c  , 


that  is , 


^(0=1  e-5t  +  |  c-2t,  Xj(t)  =  .  4  „-5t  +  4  ,-2t. 


3 •  Theory  of  a  particle  moving  from  state  to  state 

We  will  give  a  brief  development  of  the  probability  theory  that 

describes  the  behavior  of  a  particle  which  moves  from  compartment  to 

compartment  in  a  multi -compartment  system.  Wfc  assume  that  the  time  that 

the  particle  spends  in  the  ith  compartment  has  an  exponential  distribu- 

-a^x 

tion,  say  with  density  a^e  ,  and  that  the  probability  a  particle 
leaving  the  i1,1  compartment  will  go  to  the  j  compartment  is  a  constant, 
independent  of  time,  say  P^j •  Let  B(t)  =  (b^j(t))  be  the  transition 
matrix  of  the  system,  which  is  defined  as  follows:  Suppose  there  are  n 
compartments  in  the  system.  Then  B(t)  is  an  n  x  n  matrix  whose 
elements,  b^(t),  are  functions  of  t;  they  give  the  probability  that 
a  particle  placed  in  the  i  compartment  at  time  zero  will  be  in  the 
j*'*1  compartment  at  time  t.  The  key  to  the  theory  lies  in  the  fact  that 
the  matrix  function  B(t)  satisfies  B(s  +  t)  =  B(s)B(t) ,  and  hence 
there  is  a  constant  matrix  A  such  that  B(t)  =  e^.  The  argument  runs 
along  these  lines  -  to  show,  for  example,  that  B(9)  =  B(5)B(4),  and  in 
particular,  say  when  n  =  3,  that 

b23(9)  =  b21(5)b13(4)  +  b22(5)b23(4)  +  L23(5)b33(4) 

we  note  that  b23(9)  Pr°bability  that  a  particle  placed  in 


compartment  II  will  be  in  compartment  III  after  9  time  units,  say  hours. 
Now  the  only  way  that  a  particle  can  start  in  II  and  be  in  III  after  9 
hours,  is  at  t  =  5,  to  be  in  I  and  then  go  to  II  during  the  next  four 
hours,  or  at  t  =  5,  be  in  II  and  go  to  III  in  the  next  four  hours,  or 
at  t  =  5 ,  be  in  III  and  also  be  in  III  after  four  more  hours.  Thus  the 
relation  we  want  is* 

(1)  b23(9)  =  b21(5)b13(4)  +  b22(5)b23(4)  +  b23(5)b33(4)  . 

The  same  argument  applies  to  any  breakdown  of  a  time  interval  into  two 
parts,  and  hence  the  relation 

(2)  B(s  +  t)  =  B(s)B( t) , 

which  expresses  all  relations  of  the  type  (l)  simultaneously,  by  way  of 
matrix  multiplication. 

From  (2),  we  see  that  for  integers  m,  B(m)  =  B(l)m,  and  that 
B ( “)  =  B(l)^m.  Hence  for  rational  r  >  0,  B(r)  =  B(l)r  and  since 
the  elements  of  B  are  continuous  functions  of  t,  we  have  B(t)  =  B(l)^ 
for  all  t  >  0.  Writing  B(l)  =  eA,  we  get  the  basic  formula  for  the 
transition  matrix  of  the  system: 

B(t)  =  etA, 

where  A  is  a  matrix  of  constants  whose  elements  we  must  find.  We  do 
this  by  studying  the  limiting  behavior  of  B(t)  as  t  0.  Let  the 
elements  of  A  be  We  have 

(3)  B(t)  =  etA  =  I  +  tA  +••• 

*  In  this  argument,  it  is  essential  that  the  time  in  each  compartment  have 
the  exponential  distribution,  in  order  that  the  future  time  in  a  compart¬ 
ment,  given  that  a  particle  is  in  a  compartment  at  time  t,  be  independent 
of  t.  The  exponential  is  the  only  distribution  with  this  property. 


(, 


~  -a^  t  __ 

When  t  is  small,  b^  =  e  =•  1  -  a^t,  while  from  (3)  »  b^(t) 

=  ]  +  ta^  ^ .  Thus 

aii  =  "  V 

A  similar  argument  will  show  that 


aij  "  Pijai’ 


with  p .  the  probability  that  a  particle  leaving  the  ith  compartment 
jth 


A  = 


compartment.  Thus 

"al 

S1P12  a2P13 

... 

*ipA 

l2P21 

• 

"a2  a2p23 

•  •  • 

a2p2n 

• 

• 

npnl 

anpn2  anpn3 

•  •  « 

J 

We  summarize  the  theory  in  this  form: 

Theorem.  If  a  particle  moves  from  compartment  to  compartment  in  a  system 

"t»h 

of  n  compartments ,  and  if  the  time  a  particle  spends  in  the  i 

-a.  x 

compartment  has  the  exponential  density,  a^e  1  ,  and  if  p„  is  the 
probability  a  particle  leaving  the  l^1  compartment  will  enter  the 
compartment,  then  the  matrix  B(t)  ,  whose  elements  b^(t)  give  the 
probability  that  a  particle  starting  in  the  i^1  compartment  will  be  in 
the  compartment  at  time  t,  satisfies 


B(s  +  t)  =  B(s)B(t) 
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and  henoe 

B(t)  =  etA, 

where 


As  an  example,  consider  this  system. 


The  density  functions  for  the  time  a  particle  spends  in  I,  II,  and  III 
are  lle-^\  5e-"^,  and  Ve-^.  The  probability  that  a  particle  leaving 
I  will  go  to  II  is  p12  =  similarly,  p13  =  ^ ,  p21  -  tjj,  P2?  -  jjjp 


Thus 


n  -  11  -  JL 

p31  21’  p3?  ?1‘ 


/-u 

'# 

7<lf) 


11(I> 


-  5 

7(£) 


li(|>) 

.(11) 

A?.y 


and  the  transition  matrix  of  the  system  is  I3(t)  =  e 


tA 


The  above  method  may  also  be  applied  to  systems  with  absorbing 
states,  that  is,  compartments  with  flow  in  but  no  flow  out.  This  is 
done  by  putting  a^  =  0  for  such  states.  For  example,  in  this  system, 


we  put  ax  -  9,a2  =  6,&3  =  0,  p^  =  1,P13  =  O.p^  -  ^P^  =  g  and 
the  transition  matrix  is  B(t)  =  etA,  where 


/ 


A  = 


-  9 


9 

•  6 
0 


6#i 


There  remains  the  problem  of  finding  the  matrix  etA  when  the 
matrix  A  is  given.  This  may  be  done,  at  least  in  the  case  when  A 
is  similar  to  a  diagonal  matrix,  and  most  practical  problems  seem  to 
fall  in  this  category,  as  follows:  Let  the  distinct  characteristic 
roots  of  A  be  ri,r2»  *  •  *  >r,n-  Let  the  principle  idempotents  of  A 


be  E. ,E„, . . . ,E  ,  where 
1  2  m' 


Ki  = 


Jk 


(A-rU) 


9 


For  a  simple  numerical  example,  let  A  be  the  matrix  of  the  system 
given  above, 


The  rootB  are  found  by  solving 

-9A  9  0 

|A  -  X 1 1  =  2  -6 A  0  =  -K [ (9  +  X)(6  +  \)  -  18]  =  0. 

0  0  A 

The  roots  are  r^  =  -3,  r^  =  -12,  r^  =  0. 
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Hence  , 

1 1 

T 

4 

'  3 

3 

r  _  (A  +  121)  (A  -  01)  _ 

2 

2 

8 

1  (-3  +  l?K-3  -  0) 

9 

3 

~  9 

1° 

0 

0 

1 

2 

1 

3 

4 

_  (A  +  31)  (A  -  01) 

E2  (-12  +  3)  (-12  -0) 

2 

'9 

1 

3 

1 

"9 

1° 

0 

0 

\ 

0 

0 

1 

_  (A  +  31) (A  +  121)  _  1 
3  (0  +  3 )  ( 0  +  12) 

0 

0 

1 

1 

0 

0 

and 


tA  _  -3tt  -12t„  .  0tr, 

e  =  e  E^  +  e  E,,  +  e  E^ . 


If  we  write  this  as  a  single  matrix,  we  have 


B(t)  =  e 


tA 


f 


1  -3t  2  -12t  -3t  -12t 

je  +  je  e  -  e 


2  -3t  2  -12t  2  -3t  ^  1  -12t 

— e  —  — e  —e  +  Tre 


V 


o 


0 


,  4  -3t  1  -12t 

1  -  je  +  je 


,  8  -3t  1  -12t 

1  -  r  -  r 


Thus,  the  probability  that  a  particle  starting  in  II  will  be  in  I  at 

The  probability  it  will  be  in  III  Is 


+4  *  4  2  -3t 

time  t  is  ^e 


1  -r 


8  -3t  1  -12t 


2  -12t 

r  • 

If  a  trace  amount  of  substance  is  put  into  compartment 


J,  then  at  time  t,  the  proportion  of  that  amount  in  compartments  I, II, 
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and  III  will  be 


1  -3t  ,  2  -12t  -3t 

je  +  je  ,  e 


-12t 


and  1 


4e-3t  +  1^-lZt 

r  +  t 


4 .  Theory  based  on  a  particle's  sojourn  time 

In  this  section  we  will  develop  a  method  for  interpreting  the 
disappearance  of  tracer  amounts  of  a  substance  from  a  single  compartment 
or  pool.  The  method  is  different  from  those  of  Sections  2  and  3}  it  is 
more  general,  and  easier  to  apply.  The  idea  is  roughly  as  follows:  A 
radioactively  labelled  particle  enters  the  compartment  under  investi¬ 
gation.  It  spends  a  random  amount  of  time,  presumably  with  an  exponen¬ 
tial  distribution,  in  the  main  compartment.  It  then  escapes  from  the 
main  compartment,  after  which  there  are  two  possibilities  -  the  particle 
does  not  return,  or  else  it  returns  after  some  random  time  whose  distri¬ 
bution  we  wish  to  find.  We  will  call  this  the  sojourn  time  of  the  particle. 
We  consider  how  to  relate  the  tracer  disappearance  curve  to  these  two 
factors  -  the  probability  that  a  particle  will  not  return,  and  the  distri¬ 
bution  of  sojourn  time  for  those  particles  that  do  return.  Having  the 
distribution  of  sojourn  time,  we  will  consider  physical  situations  con¬ 
sistent  with  that  distribution. 

This  is  the  situation:  We  add  a  unit  amount  of  radioactive  tracer 
to  a  compartment,  and  observe  the  amount  still  in  the  compartment,  say 
f(t)  ,  at  times  t^,tg, . . .  jt^.  We  assume  that  a  radioactive  particle 
spends  a  random  amount  of  time  in  the  compartment,  and  then  escapes. 

Having  escaped,  with  a  certain  probability  it  will  not  return,  otherwise 
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it  will  return  after  a  random  sojourn  time  with  density  function  g(x) . 

The  relation  betwe*  n  f  and  g  is  given  by  this  rather  formidable 
appearing  integro-differential  equation: 

(4)  =  -  hf(t)  +  h  J  f(x)g(t  -  x)dx  . 

0 

We  will  go  through  an  example  to  show  why  (4)  describes  the  situation  and 
to  show  that  the  analysis  of  (4)  is  not  as  difficult  as  it  first  appears. 

Suppose  that  we  have  a  pool  which  contains  2  mg.  of  the  substance  in 
question.  We  use  a  radioactive  isotope  to  label  a  tracer  amount  of  the 
substance  and  inject  it  into  the  pool.  We  then  observe  the  disappearance 
curve  f(t),  the  (relative)  amount  of  tracer  still  present  at  time  t. 
Suppose  that  the  initial  slope  of  the  disappearance  curve  is  -  2.  This 
means  that  4  mg .  of  the  substance  are  leaving  the  pool  each  hour,  and, 
presumably,  4  mg .  are  entering  the  pool  each  hour,  on  the  assumption  that 
we  have  a  steady  state  : 


4mg/h 


4mg/hi 


Wow  if  none  of  the  radioactive  particles  were  fed  back  to  the  main  pool, 
the  disappearance  curve  would  satisfy 

=  -  2f  ( t) 
at 
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— 2t» 

and  hence  f(t)  =  e“  ,  since  f(0)  =1.  But  suppose  that  there  is  a 
feedback  mechanism  which  we  characterize  by  a  function  g(x).  Then  (4) 
becomes 

(5)  ~  2f(t)  +  2  5  f(t)g(t  -  x)dx, 

dt  0 

the  first  term  on  the  right  indicating  that  tracer  is  lost  from  the  pool 
at  the  rate  of  two  times  the  amount  present;  the  second,  integral  term, 
reflecting  the  cumulative  feedback.  To  illustrate  the  feedback  term, 
say  when  t  =  8  hours,  we  have 
8 

5  2f(x)g(t  -  x)dx  =  2f(l)g(7)  +  2f(2)g(6)  +  2f(3)g(5)  +•••+  ?f(7)g(l) , 

0 

the  terms  on  the  right  showing  that  at  t  =  8  hours,  the  feedback  is 
made  up  of  the  amount  received  seven  hours  ago,  2f(l),  times  the  pro¬ 
portion  of  that  amount  which  is  to  be  returned,  g(7) ,  plus  the  amount 
received  six  hours  ago,  2f(2),  times  the  proportion  of  that  amount 
which  is  to  be  returned,  g(6) ,  etc. 

Next,  suppose  that  we  have  observed  the  disappearance  curve  f(t) 
at  enough  time  points  to  have  determined  that  it  has  a  particular  functional 
form,  say 

f(t)  =  .6e"3t  +  .4e_,5t. 

The  initial  slope  of  f  is  still  -  2,  but  it  became  evident  after  a 
short  time  time  that  the  disappearance  curve  was  not  a  single  exponential, 
and  hence  that  some  of  the  tracer  was  being  fed  back.  To  find  the  distri¬ 
bution  of  sojourn  time  for  particles  returning  to  the  compartment,  we 
must  solve  (5)  for  g(x).  Now  a  little  fiddling  with  equation  (4)  will 


1/. 


establish  the  following  rule!  If  f  is  a  linear  combination  of 
exponentials ,  then  g  is  a  linear  combination  of  one  fewer  exponentials. 


In  our  case,  we  assume  that  g  has  the  form  g(x)  =  ce 


-bx 


and 


try  to  solve  (5).  The  procedure  is  elementary,  we  merely  substitute 
-bx 


ce 


for  g  and  get  conditions  on  c  and  b.  We  have 


-  1.8e~3t  -  .2e',5t  =  -  1.2e_3t  -  .8e“,5t  +  25  ( .6e"3x  +  .4e”*5x) ce“b( 1  x)dx. 


Integrating  and  collecting  terms,  this  reduces  to 


-  .3e-3t  +  .3«-5t  = 


-btv 
e  )  + 


b 


^_(e",5t 


We  thus  need  to  have 


and 


or  b  =  1.5, 


.3,  or  c  =  .75. 


Thus  our  feedback  mechanism  is  described  by 
g(x)  =  .75e~1,3x,  x  >  0. 

Notice  that  g  is  not  a  proper  probability  density  function, 
the  area  under  g  is  .5,  not  1: 

CO  oo 

5  g(x)dx  =  5  • 

0  0 


75e"1-5xdx  =  .75/1.5  =  .5. 


since 


This  means  that  only  one-half  of  the  particles  escaping  from  the  main 
compartment  are  being  fed  back}  the  other  half  do  not  return.  Thus  we 
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have  the  stochastic  interpretation  of  this  particular  disappearance 

curve,  f(  t)  =  .6e  4  '  :  A  radioactive  particle  entering  the 

main  compartment  spends  a  random  time  in  that  comnartment,  the  distrlbu- 

-2t 

tion  of  that  time  is  exponential  with  density  function  2e  .  After  a 

particle  escapes  from  the  main  compartment,  with  probability  one-half 

it  does  not  return,  and  with  probability  one-half  it  returns  after  a 

“1  •  5  x 

sojourn  time  with  exponential  density  1.5a 


We  now  look  for  physical  situations  which  can  produce  this  sort  of 

feedback)  fifty  percent  not  returned,  fifty  percent  returned  after  a 

-1  5x 

random  delay  with  density  1.5e  .  Unfortunately,  and  this  is  typical 

of  actual  problems  of  this  type,  there  are  very  many  plausible  situations 
which  will  produce  the  required  feedback  and  hence  the  given  tracer 
disappearance  curve.  For  example,  this  two-pool  system  would  have 

4-  __  CX 

tracer  disappearance  curve  f(t)  =  ,6e-  +  «4e 


But  this  pool  system  would  also  have  the  same  disappearance  curve, 


i  n 


2  mg 

4  mg/ 

2§mg 

2  mg/hr 

2  mg/hr 


2  mg/hr 


for  it  meets  the  feedback  requirements  -  half  of  the  particles  leaving 
do  not  return  and  those  that  do  return  have  a  sojourn  time  with  distri¬ 
bution  1 . ' JX. 

We  can  assign,  in  fact,  an  infinite  number  of  two-compartment 

-3t  -,5t 

systems  which  will  have  disappearance  curve  ,6e  +  ,4e  *  .  For 

any  assignment  of  r,  0  <  r  <  2,  this  two-compartment  system  will  have 
the  required  disappearance  curve: 


n 


2+r 

f(2+r) 

2 

2 

,2—r 


For  if  a  particle  leaves  the  main  compartment,  I,  then  with  probability 
2  +  r 

— - -  it  will  go  to  the  second  compartment,  and  once  there  it  will 

4 

2  _ _ _  2  +  r  _  2 


return  with  probability 


The  product  of  these  is 


2  +  r'  t'*  ‘  "  4  2  +  r 

Thus  the  probability  that  a  particle  will  be  fed  back  is  one-half,  and  the 

-1 . 5x 

random  sojourn  time  has  the  exponential  distribution  1.5e  '  ,  since 

the  time  a  particle  spends  in  II  has  that  distribution. 


5.  Relations  between  the  disappearance  curve  f,  and  the  sojourn  time, 
or  feedback  function,  g 


We  have  seen  in  the  previous  section  that  if  f(t)  is  the  amount 
of  tracer  present  in  the  main  compartment  at  time  t,  then  f  satisfies 


17 


this  dlf ferentio-integral  equations 

-  hf( t)  -f  h5  f(x)g(t  -  x)dx, 
dt  Q 

where  g(x)  is  the  feedback  function  which  describes  the  probability 
distribution  of  the  time  that  particles  spend  after  leaving  and  before 
returning  to  the  main  compartment.  The  area  under  g(x)  gives  the 
probability  that  a  particle  will  be  fed  back,  and  when  g(x)  is  normalized 
by  dividing  by  that  probability,  it  gives  the  density  function  of  the 
sojourn  time. 

In  multiple  pool  systems  with  constant  pool  sizes  and  flows,  the 
disappearance  curve  and  feedback  function  will  be  linear  combinations 
of  exponential  functions.  We  summarize  here  the  formulas  which  give  f 
in  terms  of  g  and  vice  versa  when  f  is  made  up  of  two  or  three 
exponential  components: 


When  f  and  g  are  related  by 


=  ~  hf^  +  h*  f(x)g(*  -  x)dx 

dt  0 

and  f  is  a  linear  combination  of  exponentials,  then  g  will  be  a 
linear  combination  of  one  fewer  exponentials .  Relations  between  f 
and  g  in  this  case  are  as  follows: 

1.  If  f(t)  =  p1e  1  +  p2e  2  ,  f(0)  =  l,f'(0)  =  -  h, 

then  g(x)  =  ce  ^X,  where 
b  =  a^  +  a2  -  h 

c  =  -  (h  —  a^)(h  -  a2)/h. 
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2. 


3. 


4  • 


-bx  1  ? 

If  g(x)  =  ce  ,  then  f(t)  =  p^e  +  p^e  ,  where 

and  a^  are  the  two  roots  of  the  quadratic  equation 

x2  -  (h  +  b) x  +  h(b  -  c)  =0 

and  PjiPg  are  the  solutions  to  the  linear  equations 

P1  +  p2  =  1 
pla2  +  p2al  =  b- 

An  explicit  formula  for  p^  Is  p^  =  (h  -  a?)/(a^  -  a2) . 

-a.  t  -a„t  -a_t 

If  f(t)  =  p1e  +  p2e  *  +  p  e  ,  f  ( 0)  =  l,f’(0)  =  -  h, 

-b^x  -b2x 

then  g(x)  =  c^e  +  c2e  ,  where  b-^  and  b2  are  the  roots 

of  the  quadratic  equation 
2 

X  -  (a-j^  +  a2  +  -  h)x  +  P1a2a3  +  p2al&3  +  p3&la2  =  0 


and  c2‘c2  a&tisfy  the  two  linear  equations 

ala2  +  ala3  +  a2a3  -  blb2 


C1  +  C2  =  bl  +  b2  ' 


b_c,  +  bnc0  =  b,  b„ 


h 


ala2a3 


2  1  12  12  n 

-b. x  -b„x 


-a,  t  -a_t 


If  g(x)  =  c^e  +  c2e  ,  then  f(t)  =  p^e  +  p2e  +  p^e 
where  a^.a^a^  are  the  roots  of  the  cubic  equation 


-a3t 


x?  -  (b^^+bg+hjx2  +  tb1b2+h(b1+b2-c1-c2)  ]x  -  hib-jbg-Cjbp-^hj)  =  0 


and  p^,p2,p3  are  solutions  to  the  three  linear  equations 

pl  +  p2  +  p3  =  1 

plal  +  p2a2  +  p3a3  =  h 
pla2a3  +  P2ala3  +  P3Sla2  =  blV 
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6.  Quick  and  easy  solutions  to  two  and  three-compartment  systems 

In  thi3  section  we  will  show,  by  means  of  examples,  how  easy  it 
is  to  determine  the  sojourn  time  of  a  particle  which  goes  from  compart¬ 
ment  to  compartment  in  a  multiple  pool  system.  We  need  this  basic  rule: 


If  a  particle  enters  a  compartment  which  has  several  possible 
exits,  such  as  this: 


(a+b+c)/hr 


V  units 


J 

c/hr 


a/hr 

b/hr 


a 

then  the  probability  it  will  leave  by  the  "a"  exit  is  ,  by 

the  "b"  exit,  —  b— ■  ,  etc.  Given  that  the  particle  leaves  by  the 
"a"  exit,  the  time  the  particle  spends  in  the  compartment  has  the 
exponential  distribution  with  mean  and  density  function 

e-x^a+b+c^/V.  Particles  leaving  by  the  "b"  or  "c"  exit  have 
the  same  distribution  time  in  the  compartment,  that  is,  the  distribution 
of  the  time  that  a  particle  spends  in  a  compartment,  given  that  it  leaves 
by  a  particular  exit,  does  not  depend  on  the  particular  exit,  but  only 
on  the  volume  and  total  flow  out  of  the  compartment. 


For  example,  in  this  system 


2°  T 

6 

4 

4 

lio 
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the  probabilities  a  particle  will  leave  by  the  "6",  "4",  or  "10"  exits, 
are  and  respectively,  and  particles  leaving  by  any  of 

these  exits  will  have  remained  in  the  compartment  a  random  time  with 
mean  j  and  density  function  5e  ^X,  since  =  5* 

As  an  illustration  of  how  to  find  tracer  disappearance  curves  and 
sojourn  times  for  a  system,  consider  this  two-pool  system: 


I 


n 


14 


10 

5 

9 

4 

Mm 

6 

lio 

We  ask  such  questions  as:  What  is  the  sojourn  time  of  particles  leaving 
17  What  proportion  of  the  particles  leaving  I  never  return?  What  is 
the  tracer  disappearance  curve  from  17 

First  note  that  the  probability  a  particle  leaving  I  will  go  to  II 

is  =  .5.  Once  in  II,  the  particle  will  be  fed  back  to  I  with 

probability  £—•  =  .6.  Hence  the  probability  that  a  particle  leaving  I 

will  be  fed  back  is  the  product  of  these,  or  (.5)(*6)  =  .3*  Now  for 

_5x 

particles  fed  back,  the  distribution  of  sojourn  time  is  5e  ,  since 
flow  out  of  II  is  five  times  the  amount  present.  Thus  the  feedback 
mechanism  is  described  by 


g(x)  =  ce~bx  =  .3(5e  ^x)  =  1.5e  ^x 
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Referring  to  Rule  2  In  Section  5,  we  know  that  the  disappearance  curve 
will  be 


-a  t  -at 

f(t)  =  p^  +  p2e 

where  a-^>a2  are  the  roots  of 

x2  -  (4  +  5)x  +  4(5  -1.5)  =  x2  -  9x  +  14  =  0 

that  is , 

al  =  7>  a2  =  2* 

Then 


and  hence 

f(t)  =  .4e"7t  +  .6e-2t 


is  the  disappearance  curve  of  a  unit  amount  of  tracer  placed  in  compartment 
I  at  t  =  0. 


Now  consider  another  type  of  problem:  Given  that  the  tracer 
disappearance  curve  from  a  compartment  is  a  particular  function,  say 

f(t)  =  |e"3t  +  |e'7t,  f'(0)  =  -  6, 

construct  pool  systems  with  that  disappearance  function. 

To  do  this,  we  find  the  sojourn  function,  g(x)  =  ce  ,  where, 
according  to  Rule  1  of  Section  5,  b=3+7-6=4,  and 

.  -  -  Thus 

g(x)  =  |e_6x  =  g<4e-^)  ,  X  >  0. 


1  V 

Thus  the  area  under  g(x)  is  g-,  and  hence  g  of  the  particles  leaving 
the  compartment  do  not  return;  of  the  g-  that  do  return,  the  sojourn 
time  has  an  exponential  density,  4e-4  .  There  are,  as  usual,  an 
infinite  number  of  two-compartment  systems  with  the  given  disappearance 
curve  and  sojourn  function.  For  any  choice  of  r  in  0  <  r  <  5.25, 
this  system  will  be  such  a  one: 


I 

n 

(6-r)  V 

5.25V 

W 

(M  v 

(5.25- r)Vr 

Y 

V  4  I  V 

'  .75  V 

1 

rV 


